Stability of smectic phases in hard-rod mixtures 
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Using density-functional theory, we have analyzed the phase behavior of binary mixtures of hard 
rods of different lengths and diameters. Previous studies have shown a strong tendency of smectic 
phases of these mixtures to segregate and, in some circumstances, to form microsegregated phases. 
Our focus in the present work is on the formation of columnar phases which some studies, under 
some approximations, have shown to become thermodynamically stable prior to crystallization. 
Specifically we focus on the relative stability between smectic and columnar phases, a question 
not fully addressed in previous work. Our analysis is based on two complementary perspectives: 
on the one hand, an extended Onsager theory, which includes the full orientational degrees of 
freedom but with spatial and orientational correlations being treated in an approximate manner; 
on the other hand, we formulate a Zwanzig approximation of fundamental-measure theory on hard 
parallelepipeds, whereby orientations are restricted to be only along three mutually orthogonal 
axes, but correlations are faithfully represented. In the latter case novel, complete phase diagrams 
containing regions of stability of liquid- crystalline phases are calculated. Our findings indicate 
that the restricted-orientation approximation enhances the stability of columnar phases so as to 
preempt smectic order completely while, in the framework of the extended Onsager model, with full 
orientational degrees of freedom taken into account, columnar phases may preempt a large region 
of smectic stability in some mixtures, but some smectic order still persists. 

PACS numbers: 64.70.Md,64.75.+g,61.20.Gy 



I. INTRODUCTION 

Recent work has shown that strong segregation effects 
arise in smectic (S) phases of hard rods of equal diam- 
eter but different lengths when the two components are 
mixed in varying proportions. In seminal work Koda and 
Kimura 1 analyzed a binary mixture of parallel hard cylin- 
ders using Onsager approximation, and applied a stabil- 
ity analysis. They found two types of S phases: one in 
which layers are identical and contain a mixture of both 
types of particles, and another, microsegregated phase, 
which consists of alternating layers of different compo- 
sitions. More recent work has focused on more general 
mixtures. For example, van Roij and Mulder— studied a 
mixture of parallel cylinders of different diameters using 
Onsager theory and a more complete bifurcation anal- 
ysis. They observed that, for length ratios higher than 
1:5, N-N segregation (where N stands for nematic phase) 
preempts the N-S and N-C phase transitions (C stands 
for columnar phase). More recent work by Cinacchi et 
alAi used an extended Onsager theory (EOT) on sys- 
tems of hard spherocylinders (HSPCs), but lifting the 
restriction on orientations, and complete phase diagrams 
were obtained containing I (isotropic), N, and S, phases. 
As was expected, a strong tendency toward segregation 



was observed, and microsegregated phases did also arise 
in the model, consisting of layers made up of the compo- 
nent in larger proportion -the shorter rods- alternating 
with weakly populated and highly interpenetrated layers 
of the longer rods. 

In the work by Cinacchi et al. the only spatially or- 
dered phase considered was the S phase. C phases are 
very common in nature, even in systems composed of 
prolate molecules. A simple, one-component model is 
known not to exhibit C ordering, which is preempted 
by the smectic phase. However, the C phase can be 
stabilized prior to crystallization in mixtures of HSPC 
particles with orientations restricted to be parallel, as 
shown by Stroobants et al£ using computer simulation. 
Bohle et al.fi using a density functional based on the 
weighted-density approximation, considered the effect of 
polydispersity on the stabilization of ordered phases in 
systems of parallel hard rods. Their conclusion is similar 
to that obtained previously by Bates and Frenkelji viz., 
increasing the degree of polydispersity stabilizes the C 
phase with respect to the S phase. An identical behavior 
was observed in mixtures with perfectly bimodal length 
distributions^ (i.e., two-component mixtures). 

It is intuitively clear that orientational fluctuations, 
not taken into account in all of these studies, will tend 
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to destroy columnar order, but the extent to which this 
statement is true in mixtures of hard rods is still un- 
known. As a consequence, a general concern is the ques- 
tion of whether phases of lower symmetry, C and crystal, 
may preempt all the structures observed in preliminary 
studies; in particular, all previous studies have shown the 
C phase to be stable for systems of perfectly parallel rods, 
preempting the occurrence of the smectic phase for suffi- 
ciently large length ratios. Also, in connection with the 
work by Cinacchi et al.j^. a number of questions remain 
unanswered, particularly in relation with the theoretical 
model (EOT), its particular implementation, and the ap- 
proximations used (parametrization of the distribution 
functions, etc.). These are the aims of the present pa- 
per. In this paper we have not considered the occurrence 
of crystal phases, largely for two reasons: one is that 
in mixtures the crystal phase will certainly tend to ap- 
pear at relatively high packing fractions, due to packing 
constraints, which will leave a large region where other 
spatially ordered phases may appear. Another reason is 
that the EOT is not well suited to study crystallization. 
Therefore we will limit ourselves to studying the insta- 
bility against general spatial fluctuations of the N phase. 

In the present paper we address the above questions 
from two perspectives. One is the calculation, within 
EOT, of the line where the nematic phase of the mix- 
ture gets unstable with respect to columnar-phase fluc- 
tuations. This line sets an approximate upper limit for 
the region where smectic phases should be stable. The 
other is the numerical minimization of a fundamental- 
measure functional as applied to systems of hard rods of 
square transversal section [hard parallelepipeds (HPs)] 
in the Zwanzig approximation, i.e., considering that the 
rods can assume only three possible orientations, namely, 
along the Cartesian axes xyz. In this way, by using a 
completely different approach which possesses some ad- 
vantages (better treatment of hard-core correlations and 
exact formulation of excluded-volume effects) and some 
disadvantages (restricted orientations), we may verify 
whether the predictions of Onsager theory are artifacts 
of the theory or of its particular implementation. We will 
see that some of the predictions of the two models agree, 
which gives us some confidence as to the realistic nature 
of the results. On the other hand, the results from the 
fundamental-measure theory (FMT) predict new struc- 
tures, not observed in the Onsager theory, which poses 
interesting new questions and opens up new avenues of 
research -of course the extent to which these structures 
are artifacts of the Zwanzig approximation or peculiari- 
ties of the particle's shape is as yet unknown. 

In Sec. II we briefly remind the reader the ingredi- 
ents of the two theoretical approaches. It is also devoted 
to present the methodology used to locate the spinodal 
line for C ordering in the context of the above density- 
functional theories. In Sec. IIIII we present a few results 
pertaining to the EOT used previously^ to analyze phase 
behavior in HSPC mixtures, summarizing relevant fea- 
tures that appear in the phase diagrams of a few selected 



mixtures. Complete phase diagrams, containing I, N, 
and S phases, are presented. Our purpose is to compare 
the corresponding results with those obtained using the 
Zwanzig approximation of FMT as applied to HP par- 
ticles and for equivalent mixtures. Density and order- 
parameter profiles are shown to better demonstrate the 
structure of the spatially ordered S phases. Also, spin- 
odal lines signaling the instability against C-type fluctu- 
ations are shown. The location of this line in connection 
with the stability of the S phases is discussed. We end 
in Sec. IIVI with a summary of the results and some fi- 
nal conclusions drawn from the comparison of the results 
obtained with the two theoretical approaches. 

II. THEORETICAL APPROACHES 

In this section we use two versions of the density- 
functional theory, a modified Onsager theory and 
fundamental-measure theory in the Zwanzig approxima- 
tion, to obtain phase diagrams of a few selected mix- 
tures of hard rods. Attention is restricted to spatially 
disordered phases [isotropic (I) and nematic (N)] and also 
to smectic (S) phases. The more ordered columnar (C) 
phase will be considered in Sec. II C. 

A. Extended Onsager theory for mixtures of hard 
spherocylinders 

We analyze the phase behavior of binary mixtures of 
hard spherocylinders (HSPs) of length L s (s — 1,2) and 
the same diameter a. The one-particle distribution func- 
tions of each component will be denoted by p s (r, tt). The 
density functional for the Helmholtz free energy T is writ- 
ten, as usual, as a sum of ideal .Fid and excess J- cx terms: 

F[{Ps}}=Fid[{ps}]+Fex[{Ps}}, (1) 

with the ideal part given by 

2 

{3J 7 id[{p s }] = J2 //drd«p s (r,«){lnp s (r,n)-l}, 

s=l •* 

(2) 

where (3 = 1/kT. The excess part is given by 

/?F ox [{p s }] = / dr$ ox (r;{p s }), (3) 
Jv 

where the excess free-energy density per unit thermal en- 
ergy, < I ) ex(r; {p s }), is written, in the EOT, as 

4> cx (r; { Ps }) = *MEE / dr ' II dtldtl ' 

s=l t=l ^ V "* "* 

x Ps (r, n)p t (r',n')f st (r - r', H, 

f s t are overlap functions for the three different interac- 
tions (unity if particles overlap and zero otherwise) , and 
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^(ry) is a prefactor that depends on the mean packing 
fraction r\ and that is chosen to make the theory recover 
the second virial coefficient exactly and approximate the 
remaining virial coefficients in terms of the second coeffi- 
cient, in the manner proposed by Parsons- and Lee£ Our 
implementation of the present theory (see details ini^) 
involves first using the decoupling approximation, 



p s (r, SI) ps p s (r)h s (Sl) 



(5) 



(i.e., particle positions and orientations, characterized, 
respectively, by the distribution functions p s and h s , are 
assumed to be decoupled at the level of the one-particle 
distribution functions), and then parametrizing the two 
distribution functions in terms of suitably normalized ex- 
ponential functions: 



p s (r, SI) ks p s (r)h s (Sl) 



x s p 



= A S cos qz 

lo(K) 



A s P 2 (cos( 



2tt 



dxe 



A s P 2 (x) 



,(6) 



where q = 2%/d is the wave number of the smectic un- 
dulation with period d, A s are variational parameters 
describing the orientational order (A s = in the I phase 
and A s 7^ in the N and S phases), {X s } are variational 
parameters describing the positional order (X s — in the 
N phase), x s the molar fractions of species s, p the mean 
density, and Iq(x) is the modified Bessel function of or- 
der 0. Equation (JTJ should be minimized with respect 
to the variational parameters {A s }, {A s }, and d to find 
the equilibrium density profiles and the free energy for 
particular values of p and x. The Gibbs free energy per 
particle g is obtained from 
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dp 



1 



F 



(7) 



where P is the pressure and F is the equilibrium 
Helmholtz free energy, or directly from an appropriate 
derivate of the free-energy density functional evaluated 
at the equilibrium profiles. Applying a Maxwell double- 
tangent construction on g(x) for fixed P allows us to find 
the coexistence values of the different phases in the usual 
way and to obtain the phase diagrams in the pressure- 
composition plane. 



B. Fundamental- measure theory for mixtures of 
hard parallelepipeds in the Zwanzig approximation 

The FMT used here has been described in detail in 
Ref. 10, so that only a brief sketch is provided here in 
order to show how it can be numerically implemented 
to calculate phase diagrams of mixtures. We consider a 
system of HPs, i.e., hard rods of rectangular shape. We 



continue to use the notation introduced in Sec. II A for 
the dimensions of the particles; in this case we denote 
by L s the length of a parallelepiped of species s, with 
a s being the breadth of the particles (assumed to have a 
square section). 

We have used the fundamental-measure functional 
for HP in the approximation that orientations are re- 
stricted to lie along three mutually orthogonal direc- 
tions (Zwanzig approximation), taken to be the xyz axes. 
Since the unit vector SI only has three possible orienta- 
tions, one can think of the one-particle distribution func- 
tions p s (r, SI) as being expressed in terms of a set of 
Dirac-delta functions along the three axes: 



p s (r, SI) = Ps^(r)S(Sl 



(8) 



where e M , p = 1,2,3, are unit vectors along the three 
perpendicular directions xyz, respectively. Note that this 
expression does not assume a decoupling approximation 
since it implies a strong mixing of spatial and orienta- 
tional degrees of freedom. The coefficients of these expan- 
sions, p sl j,(r), are the local density of species s parallel to 
the //-axe and represent quantitatively the orientational 
order in the system in the Zwanzig approximation. The 
free-energy functional will be given by Eqs. CJ-JSJ an d 
iJSJ, where the excess part of the free-energy density in 
reduced thermal units, obtained in Ref. 10, has the form 



$ ex (r; {p s }) = -no ln(l - n 3 ) 



ni ■ n 2 , n 2x n 2 yn2z 



1 - "3 (1 - n 3 ) 2 ' 



(9) 



with the functions {n a } (a = {0, lx, ly, lz, 2x, 2y, 2z, 3}) 
being weighted densities obtained as 



n a (r) 



where u4^ are characteristic functions whose spatial in- 
tegrals give the fundamental measures of the particles 
(edge length, surface, and volume). In this case, instead 
of the parametrization Eq. ©, used in the EOT, we use 
a Fourier expansion to represent smectic order: 



EE 

5 — 1 f-t— 1 



dv'p slx {v')^ 



(10) 



p s/J (r) = x s p 



' K 

E a< $ cos qkz 

.k=0 



(11) 



(this expression is obviously also valid for the N phase), 



where {ai k J} are Fourier amplitudes (we impose = 1 

(k) 

and uniaxial symmetry, reflected in the condition a sx = 
aiy 1 ). The coefficients {js^} (which play the role of f s {Sl) 
in the Onsager theory) are the spatially average proba- 
bilities to find a particle of species s parallel to the p 
axis (again uniaxial symmetry implies 7 SX = y s y)- The 
number of terms in the expansion, K + 1, is chosen so 

that aifP < 10~ 7 . 



(0) 
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In the nematic limit p s (r,fl) — p s h s (Cl), and the co- 
efficients 7 SM are easily obtained in terms of the nematic 
order parameters {Q s } (with —1/2 < Q s < 1) by impos- 
ing the conditions 

1 = J dflh s (fl), Q s = J dClh s (n)P 2 {n ■ e z ), (12) 

is the second-order Legendre polynomial] which 

lead to 



7sAl = -[l + (3^ z -l)Q s ] 



(13) 



with 8fj, v the Kronecker delta. We have selected the ne- 
matic director to be parallel to the z axis. 

With these definitions the weighted densities can be 
easily calculated from Eqs. ((TUll and ltTT|) as 



and the phase diagram is obtained as indicated in Sec. 
Ill Al The results for phase diagrams of different mixtures 
(including the one-component limit) pertaining to this 
theory, and also to the EOT approach outlined in Sec. II 
B, will be presented in Sec. IIIII Also, density and order- 
parameter profiles for the smectic phases will be shown 
in order to better grasp the structure of these phases. 
These profiles are defined as follows. The total density 
profile is 



(19) 



while the nematic order-parameter profile is defined by 



V ; 2 



Pa{z) 



(20) 



K 

six k=0 

where the functions G)i$ are the Fourier transforms of 
effective one-dimensional weights resulting from partial 
integration over x and y: 



>i o >*0=io(§K*)> 



(<2*) = <riv3s VB 7; k(J 



v, 
a 



q 



u>gHqk)=v sJl (1^)' 



(15) 



with v s = L s a 2 the particle volumes, cr*„ = <j s + 
(L s - er s )<5 M „, j (x) = cosx, and j 1 (x) = sinx/x. 

The ideal part of the free energy density in reduced 
thermal units for this model is 



$id(z) = 'J2ps I x(z) [lnp slx (z) - 1] . 



(16) 



Thus the total free energy per unit volume can be calcu- 
lated as 



(3T _ 1 
~V~ ~ d 



[ dz[$ id (z) + $ ox (z)] 
Jo 



(17) 



Now Eq. I|17(l is minimized with respect to the parame- 
ters 7 S(L1 , ai k J , and d to find the equilibrium density pro- 
files for fixed p and x. The pressure can be evaluated 
from 



l-n 3 (*) [l-n 3 (z)Y 



+ 2 



n2x(z)n 2 y{z)n 2 z(z) 
[l-n 3 («)] 3 



(18) 



C. Instability against spatial fluctuations 

In this section we present the formalism used to ana- 
lyze the instabilities against spatial fluctuations. We first 
describe the bifurcation analysis that we have used to lo- 
cate the instability points of the N phase against spatial 
ordering. Our aim is to see whether the phase behav- 
ior of the above HSPC mixtures is strongly modified by 
the presence of spatially ordered phases with a symmetry 
lower than that of the S phase. In the present paper we 
are interested in locating the spinodal line corresponding 
to C ordering, since this is the phase that, following pre- 
vious work on parallel hard rods, might become stable 
prior to crystallization. 

Instabilities against both S and C fluctuations can be 
explored at the same time by examining the response 
function of the N phase against general spatial fluctu- 
ations of a wave vector q. The free-energy change SF 
associated with fluctuations Sp s in the one-particle dis- 
tribution functions of the two components of the mixture 
can be expressed, up to second order in the fluctuations, 
in terms of the second functional derivative of the free- 
energy functional evaluated at the nematic phase: 

2 2 



Jp,(r,!!)J ft (r',n'), 



8=1 t = l 

5 2 T 



5p s (r,n)6p t (r',n') 



(21) 



The instability is signaled by the equation SF — 0. Two 
different approaches can be followed here. In one, we 
assume that fluctuations are realized at constant mean 
density p, so that the minimization of T is performed via 
a Lagrange multiplier or, equivalcntly, by choosing fluc- 
tuations Sp s that conserve the mean density. This is the 
approach to be followed in the canonical ensemble where 
the composition x and the mean density p are fixed. In 
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the other approach, general fluctuations, not necessarily 
mean density-conserving, are considered, and the second 
functional derivative is most conveniently expressed in 
terms of the direct correlation function C st evaluated at 
the nematic phase, 

s 2 m 



5p s (r,fl)dp t (r',fl') 



6{r - r')S(n - ft')5 st 



- C st (r-r',fl,fl'), (22) 



where C st is the direct correlation function evaluated at 
the nematic phase, 



c st (r-r',n,n') 



S 2 [fife: 



Pa (r,n)=p< N) (a 



(23) 



This route is followed in the calculations using the iso- 
baric ensemble where the quantities to be fixed are the 
composition x and the pressure P. Depending on the 
ensemble used the fluctuation has to be chosen accord- 
ingly. For the Onsager model we have chosen to use the 
canonical ensemble; in the spirit of the decoupling ap- 
proximation, we take 

5 Ps (v, ft) = e s e^pi N Hfl) = e s p s e^ T hi N \fl), (24) 

where hi N ^ is the orientational distribution function of 
the nematic phase, q an arbitrary wave vector, and e s 
an amplitude giving the strength of the fluctuation. The 
above fluctuation conserves the value of the order param- 
eters {Q s } as the mean density is fixed. The increase in 
free energy can be written in matrix form as 



V 



1 2 2 1 

9 X) X) e sPsT st etPt = -a T • T • . 



(25) 



8=1 t=l 



with a T = (ai,ct2), and a s = e s p s the components of 
a two-dimensional perturbation vector. The stability of 
the nematic phase against spatial fluctuations can be as- 
sessed by examining the signs of the eigenvectors of the 
2 x 2 — T matrix, with 



T st = T st (x;q,p) = — 

Ps 



c s M,fiA') 



£N) j.(N) 



(26) 



Here C st are Fourier transforms of the direct correlation 
function. In practice we proceed as follows. Instability is 
signaled by one of the eigenvalues of T or, equivalently, 
the determinant of T, becoming negative. Then, at fixed 
composition x, we look for the value of p for which the 
four equations 



detT = 0, 



d 



detT = 0, p = 1,2,3 



(27) 




FIG. 1: Phase diagram of a pure system of hard spherocylin- 



ders in the packing fraction (?7)-aspect ratio (k 



plane, 



as obtained from the EOT approach. Labels are N, nematic; 
S, smectic; and I, isotropic. 



are satisfied simultaneously (q^ — q x , q y , q z for p = 1,2,3, 
respectively). In the case of the HSPC model, the C 
phase consists of columns of particles arranged into a 
two-dimensional triangular lattice perpendicular to the 
nematic director, and the two-dimensional (xy plane) 
perturbing wave that makes the fluid unstable will have a 
q vector with non-zero x and y components. An instabil- 
ity associated with wave vectors q = (0, 0, q), i.e., a wave 
along the z direction (parallel to the nematic director), 
will signal instability against smectic fluctuations. For 
the system of HP the underlying lattice of the C phase is 
square, so that instability will be signaled by wave vec- 
tors q = (9,0,0) or (0,^,0). 

In the case of the FMT approximation the problem can 
be cast in matrix form by introducing the fluctuation 



Sp s (r,fl) = ps^es 



(n - e„) . 



(28) 



We remark that in this case the decoupling approxima- 
tion is not invoked, so the coupling between positional 
and orientational degrees of freedom, assumed in the im- 
plementation of the theory, is maintained. The parame- 
ters {e s/i } are taken to be arbitrary, as it corresponds to a 
general fluctuation. In particular, they are all taken to be 
independent, which is equivalent to saying that the val- 
ues of the order parameters {Q s } are allowed to vary as 
the perturbation is applied, in contrast with the method 
used in the Onsager model. This is done for the sake of 
convenience, and has no practical importance given that 
the condition for instability is searched for in the space 
(p, q), and the order parameters depend on the mean 
density p. The T matrix now has six dimensions: 



Ps\x 



(29) 



6 





FIG. 2: Phase diagram in the packing fraction (?7)-aspect 
ratio («) for a pure component system of hard parallelepipeds, 
as obtained from the FMT approach. Labels are N, nematic; 
S, smectic; OS, orientationally ordered solid; C, columnar; PS, 
plastic solid; and DS, discotic smectic. The symbols indicate 
the points corresponding to calculated phase coexistence. The 
numbers indicate that the phase transitions are of first order 
but that the scale is too small to be seen by the eye. 



with 



S 2 [PT C ; 



(30) 



being the relevant direct correlation function and 
C S p,tj/(q) its Fourier transform. Explicit expressions for 
this function will be presented in the Appendix. 



III. RESULTS 

Although the HSP and HP models are geometrically 
different (HSP particles have a semispherical cap which 
is absent in HP particles, and their sections are circular 
and square, respectively), and their theoretical treatment 
is quite different, the corresponding phase diagrams will 
be compared in this section, in the hope that we can 
extract some useful information about their actual phase 
behavior. This comparison can be made provided we use 
comparable aspect ratios k s ; we should take into account 
the relation 



HSPC 



L s +<7 



A. Phase diagrams for pure components 

Before showing the results for mixtures, and in order 
to have a feeling as to how the two theoretical models can 



FIG. 3: Phase diagram in the pressure P— composition x 
plane for a mixture of HSPC with the same diameter a and 
aspect ratios K\ = 4.5 and K2 = 8.0, as obtained from the 
EOT approach. Reduced units are used for length (a) and 
energy (kT, the thermal energy). The continuous lines in- 
dicate first-order phase transitions. The shaded regions are 
the two-phase regions of phase coexistence. The regions of 
stability are labeled by S (standard smectic formed by lay- 
ers identical in composition), N (nematic), I (isotropic), and 
S2 (microsegregated smectic phase with long particles located 
in the interlayer space). The dotted line is the spinodal line 
corresponding to the instability of the nematic phase with 
respect to columnar-type fluctuations. 



be compared and how definite predictions can be drawn 
from this comparison, it is convenient to recall the results 
obtained for the corresponding one-component systems 
using the same theoretical models that will be used later 
for mixtures. The phase diagrams for the pure systems 
are shown in Figs, ^and^ both in the packing fraction- 
aspect ratio plane (we should remind the definitions of 
the corresponding aspect ratios, k hspc and k hp , respec- 
tively). A quick comparison of the two phase diagrams 
indicates that the IN phase transition is qualitatively sim- 
ilar: it is of first order and, for decreasing aspect ratio, 
ends at some value of n, below which the N phase ceases 
to be stable and is superseded by the spatially ordered 
S phase or other more ordered phases. For HP particles 
the phase diagram has been analyzed in some detail using 
the FMT approach!! (see Fig. |2J; it is relatively com- 
plex, with the C phase becoming stable in a narrow range 
of aspect ratios, a discotic smectic (DS) phase and also a 
plastic solid (PS). At high packing fractions and k hp > 5 
the stable phase is an orientationally ordered solid (OS). 
The prediction of a DS phase has been confirmed by simu- 
lations of this particular models The phase diagram for 
HSPC obtained with the EOT has not been analyzed in 
such great detail (see Fig. It contains an IN transi- 
tion ending at some value of k which is very close to that 
predicted by the FMT approach for HP particles. The 
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X 

FIG. 4: Gibbs free energy per unit volume and unit ther- 
mal energy /3g as a function of composition x for the nematic 
phase (continuous line) and the smectic phase (dashed line) 
for a mixture of HSPC with aspect ratios (Li + a) /a = 4.5 
and (1/2 + u)/a — 8.0. A linear term 34.35 — 12.52x has been 
subtracted so as to better see the curvature of the curves. 
The circle is the bifurcation point where the nematic phase 
becomes metastable with respect to columnar-type fluctua- 
tions. The dotted lines are two possible free-energy branches 
for the columnar phase that could result from a full minimiza- 
tion of the free-energy functional. 

NS transition shows a negative slope, compared with the 
positive slope obtained in the FMT (at least in the range 
5 < k < 10), and it is of first order up to some value of k, 
above which it becomes second order. By contrast, the 
NS transition in the FMT approach is always of second 
order. Since neither the crystalline nor the plastic-solid 
phases have been explored in the EOT model (mainly due 
to the fact that this model is not expected, by construc- 
tion, to provide good results for phases with full spatial 
order) the comparison of the phase diagrams cannot be 
extended beyond the liquid-crystalline phases. This is 
not a limitation for our present purposes, since we would 
like to focus our investigation on these phases. However, 
simulations on HSPC indicate that at high packing frac- 
tions there appears a solid phase with full order, which 
becomes a plastic solid when the aspect ratio is low. A 
more elaborate theory (such as that due to Somoza and 
Tarazonai^) should be used to address this question. 

B. Phase diagrams for mixtures 

We have first investigated the phase behavior of a mix- 
ture of particles with (ki,K2) = (4.5,8.0) and <72/cri = 1. 
The results are contained in Figs. EHUl Throughout this 
section phase diagrams will be presented in the pressure 
p-composition x plane (by convention, we take the com- 
position of the mixture to be given by the variable x = x\, 
where species 1 is chosen to correspond to the shortest 
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FIG. 5: Phase diagram of the HP mixture in the FMT ap- 
proach, represented in the pressure-composition plane (re- 
duced units are used for the pressure). All particles have 
the same cross section a 2 , and its aspect ratios are ki — 4.5 
and K2 = 8. The continuous lines indicate the first-order 
phase transitions. The shaded areas are the regions of two- 
phase coexistence. The dashed lines indicate the continuous 
transitions between the nematic (N) and smectic (S) phases, 
or between different S phases. The dotted line is the spin- 
odal line signaling instability of the N phase against spatial 
fluctuations (abrupt change in slope is associated with rapid 
variation of nematic order parameter). S phase is a standard 
smectic phase. S2 phase consists of alternating layers of dif- 
ferent composition. DS is a discotic smectic phase. 

particles) . Figure [3] shows the results as obtained from 
the Onsager theory (this phase diagram has already been 
presented in Ref. 4). The most representative feature 
of this phase diagram is the strong segregation of the S 
phase, which largely preempts the FN transition. Two 
S phases appear: the standard S phase, with identical 
smectic layers, and the S2 phase, where layers of dif- 
ferent compositions alternate: this is a microsegregated 
phase. This phase is characterized by the two variational 
parameters {X s } having different signs, which implies the 
density distributions of the two species being shifted one 
with respect to the other by half a smectic period. Be- 
low a reduced pressure f3Pa 3 = 2.9 there is direct coexis- 
tence between the S and I phases, which ends in a S-N-I 
triple point. Above this pressure the S phase coexists 
with the S2 phase. The Zwanzig model for a comparable 
mixture is shown in Figs. [S]and[|)l Note that the pres- 
sure scale is very much reduced with respect to that from 
the Onsager theory, so that a good reference to compare 
the two phase diagrams is the location of the FN transi- 
tion, which is qualitatively similar to the previous case. 
In the Zwanzig model S segregation appears at a much 
higher pressure, and does not preempt the I-N transition. 
There is no direct S-I transition, which is superseded by 
continuous S-N transitions, and at high pressure there 
appears a first-order transition between the standard S 
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FIG. 6: A zoom of a particular region of the phase diagram 
represented in Fig. [^around the DS phase. For a key to lines 
and labels see caption of Fig. 
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FIG. 7: Density profiles of short (solid line) and long (dashed 
line) particles of coexisting (a) S and (b) S2 phases at pressure 
fUPa 3 = 0.5. The values of coexistence parameters are x = 
0.276, tj = 0.4768, and d/L 2 = 1.1642 for S phase, and x = 
0.693, 77 = 0.4185, and d/Li = 1.1498 for S 2 phase. 



phase and the two-layer smectic phase S2, the latter ap- 
pearing when the short-particle component is more abun- 
dant; this is in agreement with the predictions based on 
the Onsager model. However, an important difference 
is that the S and S2 phases undergo continuous phase 



transitions to the N phase, a feature that is absent in 
the Onsager model. This is to be expected since in the 
one-component HP fluid the N-S transition is always of 
second order. Also, in a model with restricted orienta- 
tions the S phase is largely destabilized with respect to 
the N phase, so that the region where the N phase is 
stable is considerably larger in the Zwanzig model than 
predicted by the EOT. In addition, the FMT approach 
predicts also a continuous transition between the S and 
S2 phases in the region where most of the particles corre- 
spond to the shorter component. This feature is absent 
in the EOT approximation. 

In order to understand the structure of the S phases 
in more detail, it is interesting to examine the density 
and order-parameter profiles along one smectic period d 
(some of these profiles were shown in Ref. 4 for the EOT 
model, so here we will only show the profiles correspond- 
ing to the FMT approach). This is done in Figs. 0OD In 
Fig. the profiles at a state point with reduced pressure 
pPa 3 = 0.5 on the S-S2 coexistence line are shown. As 
can be seen the S phase [Fig. Ufa)] is composed of iden- 
tical layers, with maxima in the density distribution of 
both species coinciding at the center of each layer. An 
interesting point to mention is that the position of the 
short particles is delocalized over the length of the long 
particles as can be seen from Fig. Ufa). By contrast, 
in the S2 phase, these maxima are shifted by half a pe- 
riod (an amount equal to d/2), so that the layers with 
different compositions alternate. Figure |SJ shows profiles 
at some particular state point on the S-S2 second-order 
transition line. In this point most of the particles in the 
system are short, and the few long particles present are 
evenly distributed among the layers and the interstitials. 
This corresponds to the second-Fourier coefficient aS of 
the density distribution of the long particles becoming 
zero at this state point. As x is reduced long particles 
tend to predominantly populate the insterstitials, defin- 
ing a bilayer structure. On the contrary, as x is increased, 
the long particles tend to arrange into the layers formed 
by the particles of the other component. This behavior 
is at variance with the predictions of the EOT approach, 
which imply that long particles added to a smectic phase 
made up of the short particles tend to populate the inter- 
stitials, regardless of their concentration. In Fig. [SI the 
density profiles of each species along a direction parallel 
to the director are also plotted with dotted lines. These 
profiles indicate that particles in the layers are mostly 
directed along the layer normal, whereas short particles 
tend to align in-plane in the interstitials, and a small 
amount of long particles also adopt this orientation. As 
can readily be seen, the order-parameter profile also sup- 
ports this conclusion. 

An interesting feature of the phase diagram in Fig. [5] 
is that, for a large concentration of short particles, there 
is a small region where a discotic-smectic (DS) phase is 
stable (see Fig. \G\ which is an enlargement of Fig. in 
the region where the DS phase appears). The DS phase 
was found in a one-component system of HP (Ref. 11) 
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FIG. 8: Density (solid line) and order-parameter (dashed line) 
profiles of the species with (a) short and (b) long particles. 
Densities are scaled with the corresponding particle volume 
v s . Density profiles of each species along a direction parallel 
to the director are plotted with dotted lines. This smectic 
phase corresponds to a state point in the phase diagram with 
(3Pa 3 = 0.442 and x = 0.9858. The smectic period in units 
of the small species length is djL\ = 1.3717, and the mean 
total packing fraction r\ — 0.3934. 



(see Fig. |21 which indicates that for a one-component 
system with k = 4.5 there appears a DS phase at pack- 
ing fractions r\ w 0.3) and confirmed by simulation^ 
The structure of the DS phase can be understood by ex- 
amining Fig. El We see that the phase consists of a suc- 
cession of identical layers but, in contrast with the usual 
smectic phases, particles in the layers are almost exclu- 
sively oriented in plane, with nematic order parameter 
large and negative. The interstitials contain a consider- 
able amount of particles, slightly oriented along the layer 
normal. The long-particle component tends to follow 
the majority component and populate the layers. The 
present result for the mixture indicates that as soon as a 
small amount of long rods is added to the pure fluid, the 
DS phase becomes very rapidly unstable. 

As an example of another mixture showing strong 
smectic segregation, but with one species composed on 
nonmesogenic particles, we discuss the cases K\ = 1, 
K2 = 8, and oijo\ = 1, i.e., mixtures of HSPC with hard 
spheres or HP with hard cubes (see Figs. 1101 and ITl"]) . 
Here again the two models give the same phase behav- 
ior: at high volume fractions of the nonmesogenic com- 
ponent we find an isotropic phase, and a relatively small 
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FIG. 9: Density (solid line) and order-parameter (dashed line) 
profiles corresponding to the species with (a) short and (b) 
long particles for the discotic-smectic phase. Densities are 
scaled with the corresponding particle volume v s . The pro- 
jection of the density distribution of each species along the 
direction of the director is plotted using dotted lines (note 
that these distributions have maxima at z/d — 0.5). The 
discotic smectic phase corresponds to the I-N-DS triple point 
with pPa 3 = 0.302 and x = 0.9715. The smectic period in 
units of the cross-section length is d/a = 1.2850, and the 
mean total packing fraction r\ — 0.2915. 

addition of this component to a fluid composed mainly of 
long particles destabilizes both liquid-crystalline phases. 
Again the N-S phase is of first order in the Onsager model 
but of second order in the Zwanzig theory. Apart from 
that, the two phase diagrams are qualitatively similar. 
The structure of the smectic phase is interesting, as it is 
again a clear example of a microsegregated phase: lay- 
ers rich in long particles alternate with layers made up 
of (almost exclusively) hard spheres (or hard cubes, as 
the case may be). Figures IT21 and H~3l containing density 
profiles, show this feature for the Zwanzig model (in the 
case of the EOT the segregation into well-defined layers 
is even more pronounced 4 ). Both models predict that 
the smectic phase is destabilized with respect to the ne- 
matic phase as the nonmesogenic component is added 14 
(the slope of the N-S2 phase boundary in the pressure- 
composition plane is positive). 

As a final mixture we have studied, using only 
the Zwanzig approximation, the case where the paral- 
lelepipeds and the cubes have the same volumes but dif- 
ferent breadths, i.e., different cross sections. The phase 
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FIG. 10: Phase diagram in the pressure P— composition x 
plane for a mixture of HSPC with the same diameter a and 
aspect ratios ki = 6.0 and K2 = 8.0, as obtained from the 
EOT approach. See caption of Fig. 1 for key to labels and 
lines. 




FIG. 11: Phase diagram of a binary mixture of hard paral- 
lelepipeds («i = 1) and cubes (K2 = 8) as obtained from the 
FMT approach. See caption of Fig. I^Jfor key to labels and 
lines. 



diagram corresponding to the case Li/cti = 8, L2/02 = 1, 
02/ci = 2 is shown in Fig. Due to the different 

cross-sectional areas of the particles the smectic phases 
strongly segregate at high pressures. The IN segregation 
region is also greatly enhanced. A remarkable feature of 
the phase diagram is the existence of demixing in the re- 
gion where the S2 phase is stable at low compositions; 
this phase transition ends in a lower critical point and, 
at some higher pressure, transforms into a triple point 
where the two bilayer smectic phases with moderate con- 



tent of cubes coexist with a third one which is mainly 
formed by cubes. The density profiles in Figs. Upland Hoi 
are useful to understand the structure of these phases. 
In Fig. HJj] the density distribution of the first two S2 
phases is shown. Figure ITST a - ) shows that there is a very 
strong microphase separation between the two compo- 
nents which therefore arrange in layers of almost pure 
composition. The origin of this isostructural phase sep- 
aration into two bilayer structures can be understood in 
terms of the so-called depletion effect, which has an en- 
tropic origin: the cubes added to a pure mixture of rods 
in a smectic configuration arrange themselves in the in- 
terstitials between the smectic layers, and this arrange- 
ment creates an effective attraction between the layers, 
giving rise to a first-order "condensation" phase transi- 
tion. 

The structure of the bilayer smectic phase at high pres- 
sure and with a high content of cubes is also very in- 
teresting as it illustrates another effect of the depletion 
interactions; this is shown in Fig. 1161 The thermody- 
namic state corresponds to the phase that coexists at 
the triple point, with x » 0.98. From Fig. I16f a) we can 
see that the cubes arrange themselves into well-defined 
layers; however, the rods adopt a very different configu- 
ration since they do not occupy the instertitials between 
layers in a uniform manner: an excess density of rods 
develops right at the edges of the main distribution of 
cubes. This is represented in Fig. Iltif b) by the solid line, 
which in turn has been split into parallel and perpendic- 
ular components to the director (layer normal) to show 
that the rods in the excess regions are oriented approx- 
imately with equal probability with respect to the three 
mutual perpendicular directions. We can conclude that 
the cubes act as a soft wall against which rods are piledup 
as if adsorbed. 



C. Spinodal lines for columnar order 

The discussion thus far has been focused on the phase 
diagrams containing I, N, and S phases. As discussed in 
the introduction, there are strong arguments that sug- 
gest that the columnar phase might play a role in the 
phase stability of the mixtures of hard rods. Using the 
methodology outlined in Sec. Ill CI we have calculated 
the spinodal line where the nematic phase becomes un- 
stable with respect to fluctuations involving local spatial 
order. As implicit in the presentation of the method, our 
approach does restrict the search for instability to any 
particular structure; usually however, the first instabil- 
ity (at the lowest mean density) is found for fluctuations 
against columnar-type fluctuations (this is signaled by a 
wave vector with nonvanishing x or y component), except 
in a few cases to be discussed below. 

The spinodal line is represented in the phase diagram 
of Fig. |3| corresponding to the mixture with K\ = 4.5, 
and K2 = 8, o^/ci = 1, as analyzed in the context of the 
EOT approach. As can be seen from the figure, the insta- 
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FIG. 12: (a) Density profiles of cubes (solid line) and rods 
(dashed line) for the coexisting S2 phase (/3P<r 3 = 0.37, x = 
0.543, 77 = 0.3579, and d/L 2 = 1.6422) at the I-N-S2 triple 
point shown in Fig. 1111 (b) Order-parameter profile of rods 
for the same thermodynamic state. 



FIG. 13: (a) Density profiles of cubes (solid line) and rods 
(dashed line) of the S2 phase which coexists with the I phase 
at /3Pa 3 = 0.5, x = 0.455, rj = 0.4414, and d/L 2 = 1.5139 
(see Fig. I1H . (b) Order-parameter profile of rods for the 
same thermodynamic state. 



bility occurs at relatively high pressure, which means that 
there is a relatively wide region of stability for the smec- 
tic phase and consequently, the predictions on smectic 
demixing in mixtures of HSPC are plausible. However, 
our results indicate that the presence of smectic segrega- 
tion should be taken with some caution, since it could be 
preempted by direct coexistence between the C and S2 
phases. Without a full calculation of the free energy of 
the C phase this question cannot be settled completely. 
In order to discuss this point, we plot in Fig. 0|the free- 
energy branches g(x) of the N and S phases for a mixture 
with fixed pressure f3Pcr 2 — 2.43. The bifurcation point 
is indicated by a circle. Starting at the bifurcation point 
at x = 0.1 two possible free-energy branches for the C 
phase have been represented (the C-phase branch will ex- 
ist for values x < 0.1, see Fig. |3J). One possibility is that 
the curve does not cross the smectic-phase branch: in this 
case the C phase exists as a metastable phase and the S 
phase is always more stable. The alternative situation is 
the one corresponding to the other curve, which crosses 
the smectic branch: now the C phase becomes more sta- 
ble and coexists with the S phase. This demonstrates 
that the bifurcation analysis can only provide trends as 
to the possible stabilization of the phases and that, once 
a bifurcation point is shown to exist, one cannot rule out 
any scenario. 



The FMT theory on the HP mixture shown in Fig. [5] 
predicts a spinodal line at low pressure, below the corre- 
sponding line associated with smectic ordering (there is a 
small region, for small values of x, where the smectic line 
is slightly above the columnar line which leaves a small 
region of smectic stability; see Figs. 151 and fTTfl . Clearly, 
the rich smectic phase behavior predicted by the theory 
is completely preempted by columnar ordering. At high 
composition x the spinodal line ends in the region of DS 
stability; this means that the DS phase preempts the in- 
stability of the I phase against columnar or plastic-solid 
fluctuations. 

Now we come to show the results for the mixture with 
Ki = 1, «2 = 8, and 02/eri = 1. The results from the 
EOT for the HSPC model are shown in Fig. [TT] We 
can see that the spinodal line shows the same trend: it 
increases with composition, and is located at relatively 
high pressure so that a wide region of smectic stabil- 
ity is found. By contrast, the FMT approach (Fig. \Tty 
again predicts a spinodal line against columnar-type fluc- 
tuations below the region of smectic stability, save for a 
small interval in composition at low values of x. 

Finally, we have also calculated the spinodal line 
associated with columnar fluctuations for the mixture 
with equal particle volumes but different cross sections, 
L1/&1 = 8, L2/02 = 1, and 02/01 = 2. The results 
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FIG. 14: Phase diagram of a binary mixture of cubes, Ki = 1, 
and rods, K2 = 8, with a\/a2 = 2 (i.e., with equal particle 
volumes), as obtained from the FMT approach. The pressure 
at which the three smectic phases coexist at a triple point is 
indicated by a horizontal line. See caption of Fig. [5] for key 
to labels and lines. 
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are in the same line as in the previous cases: the smec- 
tic phases are completely preempted by the instability 
against columnar fluctuations. Note that the slope of 
the spinodal line changes sign at x ~ 0.7 and then 
becomes a spinodal line for instability of the I phase 
(0.75 < x < 1) against spatial fluctuations; in this 
case the divergence of the structure factor occurs at the 
same density p, independent of the orientation of the 
wave vector, q = (0,0, q), (0, q, 0), or (q, 0, 0), which 
means that the instability could in principle be associated 
with the appearance of columnar, smectic, or solid order. 
This result is analogous to the behavior observed in one- 
component hard cubes^i for which it turns out that the 
smectic, columnar, and solid free-energy branches start 
at the same bifurcation point, the solid phase being the 
most stable phase. Only a complete density-functional 
minimization could elucidate the relative stability be- 
tween all these phases in the mixture. 



IV. CONCLUSIONS 

In this paper we have analyzed the phase behavior of 
binary mixtures of hard particles having different geome- 
tries, with a view to locating the conditions under which 
spatially ordered phases are formed. Two different theo- 
retical, somehow complementary, approaches have been 
used: one based on the standard Onsager theory for hard 
particles, the other being a fundamental-measure theory 
for hard parallelepipeds in the Zwanzig approximation, 
i.e. with restricted particle orientations. It is known 
that the tendency of the two species to macroscopically 
segregate is enhanced when the aspect ratio of the parti- 



FIG. 15: Density profiles for two of the S2 phases that coexist 
at the triple point with (3Pol w 0.388 (see Fig. E}. (a) 
Density profiles of cubes (solid line) and rods (dashed line) 
corresponding to the coexisting S2 phase with x ~ 0.148 (77 « 
0.4816, d/Li « 1.4384). (b) Density profile for the coexisting 
S 2 phase with xi « 0.0007 (77 « 0.4902, d/Li « 1.1934). 



cles is more dissimilar. For the more extreme cases this 
creates large demixing regions where smectic phases of 
different compositions coexist, even preempting part of 
the isotropic-nematic transition. Even more significant, 
in some mixtures the smectic phase exhibits microphase 
segregation, where the species demix at the level of the 
one-dimensional smectic unit cell. Microphases may un- 
dergo transitions to smectic phases with the usual homo- 
geneous layer composition. 

We have demonstrated that these mixtures tend to 
form columnar phases, even when the cross-sectional 
length of the particles is unequal, regardless of their rel- 
ative length. The extent to which columnar order is 
more or less stable than smectic order has been ana- 
lyzed using bifurcation theory. In general, the bifurca- 
tion analysis of the EOT model predicts that columnar 
order may appear; whether or not this order preempts 
direct smectic-phase segregation cannot be ascertained 
with the present methodology. However, it is clear that 
large regions of smectic stability still persist in the phase 
diagram. By contrast, the FMT approach in the Zwanzig 
approximation predicts that smectic order is almost com- 
pletely suppressed. The difference in behavior predicted 
can be traced back to the different treatments of orien- 
tational fluctuations (in the latter case orientations are 
restricted to lie along three orthogonal directions, which 
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FIG. 16: (a) Density profile of cubes corresponding to the 
coexisting S2 phase at the triple point of Fig. 1141 with com- 
position x « 0.98 (77 » 0.4860, d/L 2 ~ 1.1839). (b) Density 
profile of rods at the same thermodynamic conditions; the 
dotted line represents the density profile corresponding to the 
rods that lie perpendicular to the director, whereas the dashed 
line refers to the particles lying parallel to the director. The 
total density profile is represented by the solid line. 



APPENDIX: DIRECT CORRELATION 
FUNCTION IN THE FMT APPROACH. 

In this section we present explicit expressions for the 
direct correlation function of the mixtures. In the FMT 
approach we have, in the nematic phase, 
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The symbol * stands for convolution: f * g = 
J dr' f(r')g(r — r') and it is introduced as a shorthand 
notation. All expressions inside the angular brackets are 
understood to be symmetrized with respect to their in- 
dices s/x, tv. For example, 

(A.2) 

Also, all boldface variables are three-dimensional vectors 



enhances columnar order since the order parameter is 
usually very high). The results obtained with the two 
approaches probably set the limits of what should be ex- 
pected in the real mixtures. In this respect, given the 
lack of more reliable theoretical treatments that can give 
more quantitative answers and the tremendous difficul- 
ties involved in making theoretical progress, one should 
rely on computer simulations, an approach that is feasi- 
ble but unfortunately not pursued so far. 
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with a = 1, 2 and the products between any two of them 
are to be taken as scalar products. The expression for 
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The expressions for £0 and £3 and for the components of 
£ Q and C 2 are 
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with f3 = x, y, z. 
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The Fourier transform (7 SjLl! i„ (q) of the direct corre- 
lation function has the same explicit expression as Eq. 
(IA.1(I . except that the convolutions between different 
weights u>i" are to be substituted by the products of 
their corresponding Fourier transforms Cj^I) . These are 

^° ) (q) = n[ cos ( e ^ )/2 )" ' 

v=l 
3 

^ ) (q) = n[ 2sin (^ ) /2) /ft, 

< a) (q) = f tan (q a o$/2) £>$(q), 
^ Q) (q) = y cot ( ? „(tW/2) <&$(q), (A.6) 



where q = (91,(72,93) is the weight vector, and a^, 
a s + (L s - a s ) <5 M „. 
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